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Rapid clearance and disappearance of a tracer from the circulation challenges the 
determination of the tracer's binding potentials in brain (SPnd) by positron emission 
tomography (PET). This is the case for the analysis of the binding of radiolabeled 
[ 11 C]Pittsburgh Compound B ([ 11 C]PIB) to amyloid-p (Ap) plaques in brain of patients with 
Alzheimer's disease (AD). To resolve the issue of rapid clearance from the circulation, we 
here introduce the flow-independent Washout Allometric Reference Method (WARM) for 
the analysis of washout and binding of [ 11 C]PIB in two groups of human subjects, healthy 
aged control subjects (HC), and patients suffering from AD, and we compare the results to 
the outcome of two conventional analysis methods. We also use the rapid initial clearance 
to obtain a surrogate measure of the rate of cerebral blood flow (CBF), as well as a method 
of identifying a suitable reference region directly from the [ 11 C]PIB signal. The difference 
of average absolute CBF values between the AD and HC groups was highly significant 
(P < 0.003). The CBF measures were not significantly different between the groups 
when normalized to cerebellar gray matter flow. Thus, when flow differences confound 
conventional measures of [ 11 C]PIB binding, the separate estimates of CBF and SPnd 
provide additional information about possible AD. The results demonstrate the importance 
of data-driven estimation of CBF and SP ND , as well as reference region detection from 
the [ 11 C]PIB signal. We conclude that the WARM method yields stable measures of SPnd 
with relative ease, using only integration for noise reduction and no model regression. The 
method accounts for relative flow differences in the brain tissue and yields a calibrated 
measure of absolute CBF directly from the [ 11 C]PIB signal. Compared to conventional 
methods, WARM optimizes the Ap plaque load discrimination between patients with AD 
and healthy controls (P = 0.009). 
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INTRODUCTION 

The marker [ U C]PIB is a radiotracer of amyloid- P (Ap) plaque 
load, used clinically to diagnose Alzheimer's disease (AD) (Cohen 
et al., 2012). Despite its wide use, several challenges face the actual 
quantification of this and other Ap probes (Ducharme et al., 
2013), as recently reviewed by Kepe et al. (2013). Unfortunately, 
no dose-escalation and competitive displacement studies have 
been performed for any Ap imaging agent. Thus, the in vivo 
binding properties of these tracers have not been fully charac- 
terized (Villemagne et al, 2012), including the differences of 
binding to specific tissue types (Fodero-Tavoletti et al., 2009) 
such as estrogen-related receptors. Functionally, AD syndromes 
are associated with degeneration of specific functional networks, 
and amyloid deposition as measured with [ U C]PIB explains at 
most a small amount of the clinico-anatomic heterogeneity in 
AD (Lehmann et al., 2013). Even if [ n C]PIB binds uniquely to 
extracellular Ap plaque formations, the relevance of such binding 
may be coupled to autophagy as recent evidence suggests that Ap 
secretion and plaque formation depend on autophagy (Nilsson 
et al, 2013). Given its widespread clinical use, enigmatic binding 



properties, and the relative high cost of positron emission tomo- 
graphic procedures, better methods for the quantification of this 
tracer in relation to other clinical parameters adds to the clinical 
relevance of the tracer. 

The rapid initial clearance from the circulation suggests that 
[ U C]PIB is subject to flow-limited uptake, such that the para- 
metric analysis of the binding must take blood flow and washout 
kinetics into account. There is therefore a need to couple Ap mea- 
sures to both metabolic function and CBF deficits (Johannsen 
et al, 2000), as well as the loss of flow variability observed in AD 
(Rodeil et al., 2012). 

Three methodological problems complicate the correct quan- 
tification of binding of flow limited and rapidly metabolising 
tracers like [ n C]PIB in human brain (Moller et al, 2009). 

The first is the rapid disappearance of the tracer from the 
circulation and the consequent brief exchange with brain tis- 
sue. Conventional binding studies by equilibrium methods are 
affected by this absence of a continuing source of tracer in the 
circulation. Reference and binding regions independently clear 
tracer from the respective volumes of distribution, and the tracer 
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in the reference region no longer is a proper surrogate for the 
tracer in the circulation. The time-activity functions of differ- 
ent regions now depend on regional properties of binding, blood 
flow, and blood-brain barrier permeability, rather than on a 
common source of tracer in the circulation. 

The second is the uncertain choice of a region of reference of 
specific flow and no specific binding. The cerebellum is held to be 
little affected by amyloid deposition in AD, and the cerebellar gray 
matter often serves as reference region for measures of [ n C]PIB 
retention (Landau et al., 2013), but other regions have been con- 
sidered as well, including whole cerebellum (Joshi et al., 2012), 
and the Pons (Klunk et al., 2007; Knight et al, 2011). 

The third potential pitfall is the influence of cerebral blood 
flow and blood-brain barrier permeability differences on region- 
ally specific binding of [ n C]PIB. If the tracer clears from multiple 
compartments in a single region of interest with different quan- 
tities of exchangeable and bound tracer, it is possible that both 
flow and permeability differences can mimic or mask changes of 
binding. 

The problems facing conventional [ n C]PIB quantification are 
manifest when [ n C]PIB retention is evaluated in healthy subjects 
and patients with AD with different degrees of Af5 deposition. For 
example, it has been reported that PIB may bind differentially to 
polymorphic Af5 aggregates in some humans as well as in animals 
(Rosen et al, 2010; Ikonomovic et al., 2012). Additionally, appar- 
ent retention of PIB is evident in cerebral white matter both in 
vivo by PET (Fodero-Tavoletti et al, 2009) and in vitro by post- 
mortem auto-radiography (Klunk et al, 2004; Svedberg et al., 
2009). 

In the present study we aimed to establish a method of global 
parametric mapping of the binding potential (BPnd) of [ n C]PIB 
that would take the particular kinetic properties of [ n C]PIB 
exchange with brain tissue into account. For evaluation of these 
properties we compared tree non-invasive assays of [ n C]PIB 
binding in brain of healthy subjects and a group of patients 
suffering from AD. 

First, we applied a new model of tracer clearance, Washout 
allometric reference method (WARM), which we designed to map 
the washout of the tracer from regions of specific and non-specific 
binding. Previous approaches to the analysis of wash-out of trac- 
ers included the early "Height-Over-Area" method of Zierler 
(1965) and the recent "Hypotime" method of Moller et al. (2009), 
from both of which the current approach borrows, taking into 
account the methodological weaknesses discussed by Kanno and 
Uemura (1975). 

Second, measures of regional tracer clearances were used to 
identify a reference region of negligible specific binding. 

Third, the simplified reference region method (SRTM) 
(Lammertsma and Hume, 1996), also accounts for flow- 
dependent differences in tracer delivery, but contrary to the 
WARM method it assumes that differential equations modeling 
the radioactivity in a region of interest and a reference region 
are continuously coupled by exchange with a well defined arterial 
contribution of tracer (Lammertsma and Hume, 1996). 

Fourth, the reference region subsequently served to obtain 
parametric maps of -RPnd by means of a well established and 
clinically popular and simple ratio measure (SUVR) of the 



area under the specific retention curve (AUC), at a presumed 
optimal time range 40-60 min after i.v. administration of the 
tracer, relative to the AUC of the retention in the reference 
region. 

Altogether, we evaluated the methods that are used to estab- 
lish binding potentials, as well as the effect of flow-dependent 
correction on the binding potential of [ U C]PIB. 

METHODS 
SUBJECTS 

Six patients with AD (four women and two men) with an average 
age of 65 (SD = 7) years and moderately reduced Mini-Mental 
State Examination (MMSE) scores of 22-25 volunteered to com- 
plete the tomography. The patients were recruited from the local 
Dementia Clinic and screened by an experienced neurologist to 
fulfill the criteria for probable Alzheimer s disease. 

Eight healthy age-matched HC volunteers with a mean age of 
68 (SD = 5) recruited by public advertisement served as controls. 
They all had a normal physical and neurological examination and 
had a MMSE between 28 and 30. To exclude cognitive impair- 
ment they were furthermore examined with the Danish version 
of CAMCOG (Lolk et al, 2000). 

We obtained written informed consent from all subjects to the 
protocols approved by the Regional Science Ethics Committee 
in accordance with the Declarations of Helsinki. We previously 
reported some PET results from the same subjects (Rodell et al., 
2012;Gjedde et al., 2013). 

POSITRON EMISSION TOMOGRAPHY 

Image acquisition 

All subjects had positron emission recordings, one or two with 
[ ls O]water and one with [ n C]PIB, in the 3D mode of the ECAT 
High Resolution Research Tomograph (HRRT, CTI/Siemens, 
Knoxville, TN, USA) in a quiet room with the subjects rest- 
ing in a supine position with eyes open. One of the male AD 
patients only completed the [ n C]PIB recording. The images 
were reconstructed with 3D-OP-OSEM point spread function 
reconstruction (Varrone et al., 2009) using 10 iterations and 16 
subsets with FWHM at approximately 1.5 mm. The reconstructed 
images were corrected for random and scatter events, detector 
efficiency variations, and dead time. Tissue attenuation scans 
were performed using a rotating 68Ge source. Dynamic emis- 
sion recordings lasting 3 min (21 frames) were initiated upon 
bolus intravenous injection of [ ls O]water (500 MBq) or injection 
of [ n C]PIB (500 MBq). Catheters (Artflon and Venflon, Becton 
Dickinson, Swindon, UK) were inserted in the right radial artery 
and left cubital vein and arterial blood radioactivity was measured 
every half second for the duration of the PET scan by an auto- 
mated blood sampling system (Allogg AB, Mariefred, Sweden), 
cross-calibrated with the tomograph, and then corrected for 
external delay and dispersion. For anatomical orientation, high- 
resolution Tl -weighted MR images were obtained at 1.5 or 3T 
(GE Sigma Systems). 

Image registration and segmentation 

The summed emission recordings of [ ls O]water and [ U C]PIB 
were automatically co-registered to the individual MRI scans 
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using a six parameter affine transformation. Individual MR 
Images were co-registered to a locally generated version of 
the common stereotactic space (ICBM, Montreal Neurologic 
Institute) (Mazziotta et al., 2001) using a combination of linear 
and non-linear registrations (Collins et al, 1994; Grabner et al., 
2006). After the calculation of the final 16 mm non-linear PET- 
Talairach transformation grid, dynamic emission recordings were 
re-sampled into common coordinates. Regional -BPnd> Ri> and 
CBF measures were obtained from parametric PET image maps 
using standard model based segmentation (Collins et al., 1994; 
Grabner et al., 2006). The regions analyzed were cerebral cortex 
excluding cerebellum (CORT), putamen (PU), caudate nucleus 
(CN) frontal (FL), occipital (OL), parietal (PL), and temporal 
(TL) lobes, as well as white matter (WM), and the cerebellar gray 
matter (CERB). 

QUANTIFICATION OF [ 11 C]PIB RETENTION 
Reference region ratio measure (SUVR) 

The [ n C]PIB retention can be calculated by determining the 
accumulation relative to a reference tissue to obtain a ratio 
measure (SUVR). The ratio measure is the fraction of the region- 
of-interest integral of [ n C]PIB accumulation at steady-state, 
assumed to have been established no later than this time after 
injection (f s = 40 min), extended to the end (t e = 60 min), rel- 
ative to the integral of the [ U C]PIB accumulation observed in 
the same period in the reference region. In the reference region, 
we assume the accumulated tracer as function of time, fMND(f) 
to represent non-specific binding after delivery of the tracer by 
homogeneous flow to all voxels of the reference region. In this 
context, we further assumed the interval from t s = 40 min to 
t e = 60 min to be sufficient to establish steady-state or secular 
equilibrium in all regions, as the basis for the definition of the vol- 
ume of distribution of the tracer, Vj, as the sum of the volumes of 
distribution of non-displaceable tracer (Vnd) and an additional 
volume of distribution of displaceably bound tracer. 

Washout allometric reference method (WARM) 

In the case of negligible input from the circulation after the initial 
brief uptake, the tissue time-activity curves of the radio ligand are 
established by the radioactivity initially persed to the tissue and 
the subsequent washout from the brain regions of uptake. The 
WARM method specifically takes this condition into account and 
uses only the differences among washout rates from regions with 
different properties of binding, blood flow rates, and blood-brain 
barrier permeability. The condition means that the differential 
Equations (1) and (2) 



dm*(t) 
dt 



= K 1 c a (t)-k 2a m*(t) 



and 



df 



c a (t) - k 2 m ND (f) 



(1) 



(2) 



are linked only while the tracer is dispersed from well-defined c a , 
i.e., during the brief uptake period until maximum peak (within 
2-10 min timeframe) when washout is assumed to be negligible. 



Equations (1) and (2) also form the base of the two compart- 
ment Simplified Reference Tissue Model (SRTM) (Lammertsma 
and Hume, 1996). The term K\ is the unidirectional clearance of 
the tracer c a by the tissue, K^ D is the clearance of the tracer c a by 
the reference region, m* and mj^ D are the measured PET signal 
in the tissue (with displaceable binding) and reference, respec- 
tively. The term k 2a defines the apparent measurable washout rate 
constant for the ROI. The term k 2 is the unknown washout rate 
for non-specifically bound tracer of the same region of interest, 
and k^ D defines the measurable washout rate of non-specifically 
bound tracer in the reference tissue into the plasma. The uncou- 
pling of the first and the second term on the right hand side of the 
equations means that elimination of the first K\ and Kj terms 
yields the equations for the washout part of the signal. 



dm*(t) 
dt 

^ndW 
dt 



-k 2a m*(t) 
-k™ D m^ D (t) 



(3) 
(4) 



The rate constants k 2 and k 2 of washout of non-specifically 
bound tracer from tissue to plasma are linked by R\, 

k 2 = Rikf D (5) 

The distribution volume ratio(DVR) can be expressed as an allo- 
metric relationship between the logarithmic of the fraction of 
remaining tracer in a ROI relative to the deposited amount before 

washout, and the reference region where ^^p^ is the scaling 
exponent, found in this formulation by log-log linearization as, 



In 



ND 



(0 



'ND 



(0) 



m*(t) 
m*(0) 



(6) 



where the ratio m*(t)/m*(0) is defined by the amount of tracer 
m*(t) remaining relative to the amount of tracer m*(0) initially 
deposited before the washout. Similarly for m^ D (£)/mj^ D (0). 

For direct calculation without linearization, BPnd can be 
found using the operational equation 



BPnd(T) 



m*(0) J Q r (ln Kd(0) - In (m* ND (0))) df 
™nd(°) jo T (ln(m*(f)) -ln(m*(0))) dt 



1(7) 



Intuitively, when log transformed, the fraction part (i.e., DVR) of 
this equation states that the nominator is the accumulated log- 
signal for the reference tissue relative to how much was present 
before washout, this difference is scaled by the initial tracer 
amount of the ROI. The denominator describes the accumulated 
log-signal for a ROI or voxel relative to how much was present 
before washout. This difference is scaled by the start amount of 
the reference region. Thus the fraction is corrected both for flow, 
i.e., initially deposited tracer, and the exponential behavior of the 
washout. For a more detailed derivation please refer to the Theory 
section. 

Figure 1 illustrates the dynamics of the nominator, the denom- 
inator, and the BP^d(T) terms of Equations (43) and (40) (in the 
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Theory section) for simulated ROI and reference curves with 20% 
Gaussian noise added . As seen, the nominator and the denomi- 
nator from Equation (40) (in panel 2 from the left) are stabilized 
by the integration in Equation (43) (panel 3 from the left). Panel 
4 illustrates the convergence of the BP^u(T) estimates toward the 
theoretical result. 

Simplified reference tissue method 

The simplified reference tissue method (SRTM) (Lammertsma 
and Hume, 1996) yields binding potential when a single tis- 
sue compartment model fits the data. SRTM solves differen- 
tial equations similar to (25) and (26) (Equations (1) and (5) 
in Lammertsma and Hume (1996)). The method applies stan- 
dard non-linear regression analysis to establish the relationship 
between the tissue concentrations of the region of interest and a 
reference region tissue such as typically the cerebellar gray mat- 
ter in the case of [ U C]PIB. The parameters estimated are Rl, ki, 
and 5-Pnd- -Rl accounts for differences of delivery to the regions 
of interest and reference. 

Flow dependence of specific binding measure 

With the regional flow ratio measure R\ derived from [ ls O]water 
or [ 1 1 C] PIB analysis for each voxel with a signal n ( t ) , we simulated 
flow-adjusted reference curves (mndW) for the corresponding 
voxel of the image, determining the tracer washout by the actual 
flow measured in the voxel. Here, the term «nd(0 represents 
the simulated dynamic wash-out of the tracer that would have 
been recorded from an individual voxel in the absence of bind- 
ing. The result is a simulated image of the dynamic passages of 
the tracer through every voxel as functions of time in the absence 
of any binding in any voxel. The simulation is based on the 
dynamic behavior in a fixed reference region (cerebellar gray mat- 
ter) given by mND(0> which is the dynamic time-activity curve 
of the reference region. Figure 2 shows an example of the real 
and simulated non-displaceable (unbound) tracer time-activity 
curves for a small white matter region and a small putamen 
region, as well as the real time-activity curves for the ROIs. The 
obvious variation in the flow corrected reference curves illus- 
trate that a single uncorrected reference curve may bias the result 
significantly. Further details can be found in the Theory section. 




FIGURE 1 | Behavior of the parts of the WARM method. The first 
leftmost panel shows two initial simulated mono-exponential washout 
curves, a reference curve (k^ D = 0.92) in red, and a binding curve 
(^2a = 0.98) in blue. The black curve shows the flow normalized version 
of the reference curve normalized to the ROI curve, with 20% Gaussian 
noise added. The second panel shows the denominators and nominators 
from Equation (40) as function of time, i.e., the fraction between these 



QUANTIFICATION OF [ 15 0]Water ACCUMULATION AS CBF [k" 2 °) 

The use of positron emission tomography with [ ls O]water in sub- 
jects included in the present analysis for the purpose of correction 
of measures of CBF for the effect of CO2 tension in blood were 
reported previously (Rodell et al., 2012). We quantified the CBF 
as the unidirectional blood-brain [ ls O]water clearance (K^ 2 °) in 
units of ml hg _1 min -1 with the linearized two-compartment 
model (Blomquist, 1984) modification of Ohta et al. (1996) 
and the Lawson-Hanson non-negative least squares solution to 
general least squares functions (Lawson and Hanson, 1974). 

AUTOMATIC DETECTION OF REFERENCE REGION 

In order to validate the cerebellar gray matter as a non-specifically 
binding region, we mapped the distribution of the washout index 




two ideally linear factors form the binding potential. The third panel 
shows the similar integrated nominator and denominator from Equation 
(43). The integration stabilizes the result. The fourth panel shows the 
PBnd from Equation (40) in blue, and from Equation (43) in black. The 
dotted line is the theoretical true value for PB ND . The PB ND from both 
Equations (40) and (43) converge toward the true value with increasing 
tomographic duration. 



PIB 



14 - 




0 -1 1 1 1 i 1 1 1 1 : 1 

0 10 20 30 40 50 60 70 80 

time(min) 

FIGURE 2 | Time activity curves of PIB from a single subject in different 
regions of interest (ROI). In blue is shown the measured cerebellar 
reference curve. The purple and yellow curves are the putamen curves as 
well as the estimated reference curve. The red and green curves are the 
curves for a small white matter region and the corresponding estimated 
reference curve. 
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© image, which has unit of time, as defined for the HYPOTIME 
method (Moller et al, 2009). Distinctive regions of high 0 values 
(i.e., close to T = Tomography duration) are indicative of a poten- 
tial reference region, as previously demonstrated for the tracer 
[ n C]WAY-100635 ([ n C]WAY) (Hirvonen et al., 2007; Moller 
et al., 2009) where the cerebellar white matter served as the 
reference region. 

RESULTS 

WASH-IN AND WASH-OUT PHASES 

Due to the rapid removal of the tracer from the circulation, 
analysis of the dynamic [ U C]PIB record revealed an initial high 
frequency signal of arterial origin in the first 2 min of record- 
ing, followed first by a maximum peak (within the timeframe 
2-10 min) from which there was only by washout. For the WARM 
method, we split the signal into the three time frames of, first, 
arterial phase, second, peak uptake, and third, wash-out. We used 
the maximum peak within (the 2-10 min peak uptake phase) for 
estimation of the relative uptake coefficient P™ and the cali- 
brated surrogate CBF index directly from the [ U C]PIB signal. We 
used the wash-out phase from maximum peak within 2-60 min of 
the dynamic record for estimation of the binding potential BPnd- 

SIGNAL-TO-NOISE RATIO AND STABILITY 

In order to establish the temporal dependence and stability of the 
WARM method in relation to tomography duration, we analyzed 
the performance of the WARM and SRTM methods at different 
scan-time. After 10 min, the signal-to-noise ratio increased for 
[ n C]PIB because of the decay of the radioactivity and washout 
of the tracer. After 60 min, the standard deviation reached almost 
half the signal, as illustrated in Figure 3 (left panel). Figure 3 
(right panel) shows the stability of the SRTM and WARM method 
results with respect to time in the tomograph. The WARM 
method results converge after 60 min, while the SRTM results 
retain some dependence on time. Considering the scale used 
for BPnd> the results of the two methods are in relative good 



SNR 

4 1 i i i i i i r 




1.5 - 



0 10 20 30 40 50 60 70 80 

Time(min.) 

FIGURE 3 | The left panel shows the relative noise {SNR) of a 
typical PIB scan as a function of tomography time duration for a 
region covering the cerebral cortex without cerebellum. After 60 min 
the noise level is more than half the signal, probably due to the loss 
of the signal due to decay of ["C] and washout of the tracer, the 
right panel shows the stability of the BP ND measure for the WARM 



agreement. Considering the SNR and stability together, we regard 
a maximum time of 60 min to be an acceptable compromise. For 
this reason, we confined the calculations of binding potential with 
the WARM method to the 2-60 min time window. The SRTM 
calculations were made with a maximum duration of 60 min. 



PIB RETENTION 

In order to test the methods' power to distinguish the AD from 
the HC subject groups by quantification of the amount of reten- 
tion, we extracted regional parametric values for the PIB retention 
from the parametric images with the three different methods 
WARM, SRTM and SUVR. For the WARM, SRTM methods, we 
report mean absolute BPnd values and for the SUVR method we 
extracted DVR values. 

The mean absolute BPnd values from the WARM and SRTM 
methods were in very good agreement for cerebral cortex val- 
ues for the AD patients, while for the HC subjects the WARM 
method values are lower for most of cortex. The WARM method 
therefore yielded more significant differences between the AD and 
HC groups than the SRTM method for cortex (CORT), putamen 
(PU), frontal lobe (FL), occipital lobe (OL), temporal lobe (TL), 
and parietal lobe (PL). Only caudate nucleus (CN), white matter 
(WM) and cerebellum (CERB) reference values were less signif- 
icant than with the SRTM method. Compared with the SUVR 
method, both the WARM and SRTM methods yielded greater 
between-group differences of retention. Table 1 lists the signifi- 
cance levels, and Table 2 lists the mean parametric estimates for 
each region and group. 

Closer examination of the WM binding revealed a consider- 
able amount of white matter retention with the WARM method, 
although the flow normalization accounted for some of this. The 
SRTM method yielded negative binding values for these WM 
areas (not visible within the specified range). We judge the neg- 
ative values to be due to difficulties of regression to the WM 
dynamics. 



BP. 




20 30 40 50 60 70 



Last frame used (min) 

method and SRTM, for each method is shown the result of applying 
the analysis to images with different end frame times. The curves 
represent the parametric BPnd values for a typical AD patient for a 
single voxel in frontal lobe and a striatal region mean value. As seen 
the methods are in relatively good agreement, but the WARM method 
converge with more and more data considered. 
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Table 1 | Regional statistical significance (student t-test) between AD vs. HC subject values for both [ 11 C]PIB retention and flow estimates 
(*p<0.05) ("p< 0.001). 



Region PIB-retention Flow 



WARM SRTM SUVR R^ IB Rt CBF CBF PIB PIB "CBF" 

(H 2 0) (pC0 2 corrected) "CBF" (pC0 2 corrected) 



CORT 


♦0.0092 


•0.0206 


•0.0410 


0.0366 


0.0693 


•0.0029 


** 0.0007 


0.1933 


0.2032 


Pu 


* 0.0099 


•0.0232 


• 0.0299 


0.4069 


0.4396 


•0.0076 


• 0.0064 


0.2582 


0.3250 


CN 


•0.0161 


•0.0078 


0.4803 


0.0848 


0.1144 


•0.0088 


••0.0007 


0.0762 


0.0537 


FL 


•0.0139 


•0.0282 


* 0.0379 


0.0654 


0.2033 


•0.0045 


• 0.0016 


0.2085 


0.2343 


OL 


•0.0119 


0.0600 


0.2479 


0.0771 


0.0512 


•0.0050 


•0.0023 


0.1834 


0.2113 


TL 


•0.0026 


•0.0170 


•0.0228 


•0.0065 


• 0.0099 


•0.0016 


••0.0001 


0.1341 


0.1165 


PL 


•0.0129 


•0.0197 


• 0.0470 


•0.0191 


•0.0277 


•0.0041 


••0.0006 


0.1323 


0.1150 


WM 


•0.0241 


•0.0146 


0.7859 


0.1006 


0.1443 


* 0.0074 


•0.0011 


0.3473 


0.4690 


CERB 


0.8383 


0.3334 


0.7517 


0.3682 


0.4561 


•0.0190 


•0.0196 


0.6840 


0.8808 



Table 2 | BPnd and DVR values measured for AD vs. HC subjects using three different methods ('indicates statistical significance between 
subject groups p < 0.05). 



Region 


WARM BP 


ND 


SRTM BP ND 


SUVR DVR 


AD 


HC 


AD 


HC 


AD 


HC 


CORT 


•0.243 (0.146) 


•0.083 (0.018) 


•0.267 (0.125) 


•0.140 (0.043) 


* 1.614 (0.460) 


* 1.209 (0.165) 


Pu 


•0.280 (0.163) 


•0.101 (0.030) 


•0.422 (0.139) 


•0.270 (0.077) 


•2.025 (0.512) 


* 1.517 (0.227) 


CN 


•0.282 (0.169) 


•0.108 (0.041) 


•0.339 (0.129) 


•0.174 (0.055) 


1.519 (0.499) 


1.373 (0.229) 


FL 


•0.285 (0.192) 


•0.090 (0.023) 


•0.322 (0.175) 


•0.156 (0.058) 


* 1.790 (0.607) 


•1.248(0.210) 


OL 


•0.180 (0.097) 


•0.075 (0.021) 


0.194 (0.88) 


0.123 (0.033) 


1.385 (0.360) 


1.214 (0.147) 


TL 


* 0.204 (0.096) 


•0.074 (0.016) 


•0.222 (0.80) 


•0.129 (0.041) 


• 1.528 (0.365) 


•1.144 (0.161) 


PL 


•0.281 (0.189) 


•0.087 (0.020) 


•0.299 (0.146) 


•0.149 (0.049) 


• 1.835 (0.593) 


* 1.328 (0.217) 


WM 


•0.258 (0.073) 


•0.173 (0.045) 


•0.189 (0.065) 


•0.116 (0.026) 


1.990 (0.268) 


1.954 (0.204) 


CERB 


0.107 (0.034) 


0.111 (0.030) 


0.111 (0.009) 


0.103 (0.015) 


0.981 (0.034) 


0.987 (0.035) 



To illustrate the regional distribution of the binding, Figure 4 
shows the group mean binding images for the three different 
methods. The WARM and SRTM method results are reported as 
BPnd values, scaled similarly, while the SUVR values are reported 
as DVR values in the 1-3 range. 

CEREBRAL BLOOD FLOW 

Surprisingly, the absolute CBF values derived from [ ls O]water 
(Table 3) most significantly distinguished the AD from the HC 
groups (Table 1), especially after correction for the within-group 
variation of arterial CO2 tensions. The low standard deviations 
compared to the retention values probably accounted for the 
more significant differences among the CBF estimates. The CBF 
values of the reference region (cerebellar gray matter) differed 
among the groups. This difference disappeared when we consid- 
ered the relative flow indices, i.e., the -Ri values (Table 4). Thus, 
normalization may remove important differences and should be 
used with caution (Mayr, 1982; Borghammer et al., 2009a,b). 
The Ri estimates obtained with [ ls O]water and [ n C]PIB were 
similar (Table 4), and for both methods the temporal and pari- 
etal lobes were the only regions with significantly different values 
when we compared AD and HC subjects. Hence, it is feasible 
to use Ri estimates from the [ n C]PIB sequence to correct 



for intra-subject flow normalization, also when inter-individual 
flow differences may be of clinical value as a marker of AD 
pathology. 

We used a site-specific scale factor to obtain a simulated abso- 
lute flow estimates from the PIB peak arrival characteristics, but 
this factor was derived only from healthy subjects. We did not 
find the same statistical differences for the simulated CBF val- 
ues (Table 5) as for the standard CBF estimates. As illustrated for 
select regions in Figure 5, this failure maybe due to the high vari- 
ability observed in the simulated CBF values, as generally there 
was good correlation for the group means seen in Figure 6, except 
for white matter and cerebellar gray matter. Here, the simulated 
CBF values did not yield the same CBF value decrease for the AD 
subjects. 

REFERENCE REGION 

A clear advantage of the HYPOTIME method was the produc- 
tion of washout indices in the form of 0 images with unit of 
time as an intermediate result (Moller et al., 2009). High © 
values (0 > 1660 s relative to 1= 1890 s) are indicative of the 
region's potential for service as reference. This was previously 
demonstrated in studies of the WAY tracer, where the reference 
region was selected as the cerebellar white matter. However, for 
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Table 3 | Regional absolute CBF measures before and after correction 
for PC0 2 (Rodell et al., 2012). 



Region 




CBF 


CBFIdCO-j 


corrsctsd) 




AD 


HC 


AD 


HC 




ml hg _1 min" 


" 1 ml hg _1 min -1 


ml hg _1 min -1 


ml hg _1 min -1 


CORT 


* 40.14 (6.0) 


* 51.8 (5.0) 


** 37.6 (1.2) 


54.4 (8.0) 


Pu 


53.8 (5.6) 


65.5 (6.7) 


*50.7 (3.9) 


68.9 (11.6) 


CN 


*35.4 (6.0) 


*47.5 (7.1) 


**33.1 (1.4) 


49.6 (7.8) 


FL 


*43.7 (6.6) 


* 54.4 (4.4) 


M0.9 (2.0) 


57.3 (8.5) 


OL 


*39.0 (4.8) 


*50.3 (6.2) 


*36.7 (3.4) 


52.8 (8.6) 


TL 


* 34.3 (7.4) 


» 48.4 (4.9) 


»* 31.9 (2.0) 


50.8 (7.1) 


PL 


*39.5 (8.4) 


*54.0 (6.1) 


**36.8 (3.7) 


56.6 (8.7) 


WM 


*30.8 (4.7) 


*39.1 (4.3) 


*28.9 (1.4) 


41.1 (6.0) 


CERB 


*46.9 (4.8) 


*54.3 (4.7) 


•44.3 (4.7) 


57.2 (9.8) 



CBF values derived from the [ w 01water images. 

* and ** indicates statistical significance between subject groups (*p < 0.05) 
(**p < 0.001). 



FIGURE 4 | Images of [ 11 C]PIB binding in brain. Upper three rows show 
mean binding images of sections of the brains of a group of AD patients 
(A/ = 6), and lower three rows show similar sections for age-matched 
healthy controls (N = 8). For each group is shown the BPnd values derived 
from the WARM and SRTM methods. For the SUVR method, the mean 
DVR measures are illustrated in the range 1-3. 



the present tracer, the cerebellar gray matter consistently stood 
out as a reference region candidate. Figure 7 shows the average 
© washout index image for all subjects, for comparison with 
an insert of the reference region previously identified for WAY 

(Moller et al, 2009). 

DISCUSSION 

For [ n C]PIB, the mono-exponential washout kinetics with rapid 
disappearance of tracer from the circulation complicates the 



Table 4 | Regional /?i measures relative to the cerebellar gray matter 
reference. 



R™ 3 RtHzO 



Region 


AD 


HC 


AD 


HC 


CORT 


0.73 (0.05) 


0.82 (0.08) 


0.74 (0.06) 


0.82 (0.07) 


Pu 


0.99 (0.06) 


1.05 (0.13) 


0.99 (0.06) 


1.05 (0.14) 


CN 


0.65 (0.09) 


0.76 (0.11) 


0.65 (0.10) 


0.76 (0.11) 


FL 


0.78 (0.07) 


0.87 (0.07) 


0.81 (0.09) 


0.87 (0.07) 


OL 


0.72 (0.04) 


0.80 (0.08) 


0.72 (0.03) 


0.80 (0.08) 


TL 


•0.62 (0.08) 


*0.77 (0.07) 


*0.63 (0.09) 


•0.77 (0.07) 


PL 


*0.72 (0.09) 


*0.86 (0.09) 


*0.72 (0.10) 


*0.86 (0.09) 


WM 


0.56 (0.04) 


0.62 (0.07) 


0.57 (0.05) 


0.62 (0.07) 


CERB 


0.86 (0.2) 


0.88 (0.03) 


0.86 (0.02) 


0.88 (0.03) 



The measures were derived either from I 11 CIPIB images or from the I 15 Olwater 
images (* indicates statistical significance between subject groups p < 0.05). 



Table 5 | Regional surrogate CBF measures before and after correction 
for PC0 2 (Rodell et al., 2012). 



Region 


PIB 


'CBF" 


PIB "CBF"(pC0 2 


corrected) 


AD 


HC 


AD 


HC 


CORT 


44.8 (7.1) 


52.1 (10.2) 


44.9 (6.8) 


49.8 (6.5) 


Pu 


59.7 (7.0) 


67.7 (13.9) 


59.9 (5.7) 


64.9 (9.7) 


CN 


40.6 (6.8) 


51.2 (10.7) 


40.8 (6.1) 


48.9 (6.9) 


FL 


46.5 (7.1) 


53.3 (9.9) 


46.7 (6.3) 


51.2 (6.2) 


OL 


44.5 (7.1) 


51.8 (9.9) 


44.7 (6.3) 


49.7 (6.6) 


TL 


41.7 (7.5) 


50.2 (10.1) 


41.8 (6.3) 


48.1 (6.5) 


PL 


44.4 (9.6) 


54.0 (10.7) 


44.5 (8.3) 


51.7 (6.8) 


WM 


36.3 (4.8) 


40.6 (9.0) 


36.5 (4.3) 


39.0 (6.7) 


CERB 


57.2 (5.4) 


59.3 (10.9) 


575 (5.5) 


56.9 (6.6) 



CBF values derived from the f" CjPIB images. 
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FIGURE 5 | Distribution and mean relationships between the | 11 C| PIES derived surrogate CBF values (lower panels) and the standard [ 15 0]water 
derived CBF measures (upper panels) for selected regions in healthy controls (left panels) and AD patients (right panels). 



qualitative assessment of binding with existing methods. Here 
we demonstrate the use of the novel flow-independent WARM 
to calculate the binding potential of [ C]PIB. The calculation 
is based on the clearance of the tracer from the brain, relative 
to the initial deposits in brain tissue, compared to the same 
measure in a non-binding reference region with the same ini- 
tial deposit and washout characteristics of the cerebellar gray 
matter. We measure the binding either regionally by allometric 
linearization using Equation (41), or by direct parametric map- 
ping using Equation (43). Compared to the standard SRTM and 
SUVR methods, WARM yielded binding measures with greater 
statistical discrimination between groups of AD and HC subjects 
for cortical gray-matter (excluding cerebellum), putamen, frontal, 
occipital, temporal, and parietal lobes, and lower differences 
for caudate nucleus, white matter, and especially cerebellum, as 
expected. This has relevance for the enhanced use of quantified 
[ C]PIB retention for clinical discrimination of AD, regardless of 
the relationship between amyloid-^ deposits and the disease. 



The WARM method revealed a converging stability of the 
BPnd measure as a function of tomography duration with accept- 
able stability at 60 min, at a time when the SNR is still favorable. 
The method also takes advantage of direct calculation without 
regression. The direct calculation is important to its use because 
it serves to reduce tomography duration without compromise of 
the quantification, without the loss of the uptake signal, which we 
have shown to be useful for flow estimation. 

Based on the 0 washout index with unit of time, defined in 
Moller et al. (2009), we found that the cerebellar gray matter con- 
sistently is the best choice of reference region with no specific 
binding of [ U C]PIB. The striking difference between the cere- 
bellar reference region for tracer WAY (Hirvonen et al., 2007; 
Moller et al, 2009) and [ n C]PIB, testifies to the value of directly 
confirming the presence of a true reference area from the PET 
sequence. 

We also showed the feasibility of using the initial K\ sig- 
nal as a surrogate measure of absolute CBF, directly from the 
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Estimated vs Measured CBF 






Slope 


HC 


0.9980 i 0.09756 


AD 


1.045 1 0.1104 



40 60 
Measured (ml hg" 1 min" 1 ) 



FIGURE 6 1 Left panel shows relationship between the |"C|PIB 
estimated surrogate CBF measures and the standard 
[ I5 0]water-measured CBF values for averages of healthy controls 



and AD patients (ml hg -1 min -1 ). Right panel shows the linear 
relationship between the different regions averaged for the HC and 
AD groups. 



[ n C]PIB image. However, compared to absolute CBF measures 
with [ ls O]water , the CBF index derived from [ U C]PIB had 
greater variability. As [ U C]PIB is comparatively lipophilic, there 
are probable permeability differences from [ ls O]water images 
(Gjedde et al, 2013), which may account for some of the vari- 
ability. The absolute CBF estimates differed significantly for all 
regions, especially after correction for arterial PCO2 differences 
(Rodell et al, 2012), notably also for the cerebellar reference 
region. When we calculated the R\ values relative to the cerebellar 
reference, much of the difference was eliminated, and only pari- 
etal and temporal lobes had significant decline in the AD group. 
This observation shows that the normalization to relative flow 
measures may mask important disease specific information in 
AD, although it does enable the correction for the intra-individual 
flow bias (Borghammer et al, 2009a). The observation is of con- 
cern also to the MRI estimation of CBF indices which rarely 
are reported as calibrated measures. Clinically, the absolute CBF 
value should be regarded as an important imaging parameter for 
the diagnosis of AD, alongside the [ U C]PIB retention, the flu- 
orodeoxyglucose (FDG) uptake, and the accumulation of other 
markers. 

Non-specific binding in white matter (WM) areas is found 
both when assessed in vivo by PET (Fodero-Tavoletti et al., 2009) 
and by postmortem autoradiography (Svedberg et al, 2009). 
This implies some non-specific binding to WM, also when tis- 
sue slices are superfused with the tracer, rather than accumulating 
the tracer from the circulation after transport across the blood- 
brain barrier. However, we find that the WM binding often 
observed with [ n C]PIB is due both to slow washout kinetics and 
to non-specific retention that is not explicable by R\ differences. 

In conclusion, we demonstrate that the WARM method 
yields a stable measure of PIB's binding potential with rela- 
tive simplicity and reasonable tomography duration, employ- 
ing only integration for noise reduction with no need for 



model regression. The washout index © is readily used to 
assess the suitability of a potential reference region. The 
method accounts for relative flow differences in the tis- 
sue and yields a calibrated measure of the absolute CBF, 
obtained directly from the [ U C]PIB signal. Comparable with 
the SRTM and SUVR methods, the WARM method provides 
better discrimination between AD subjects and healthy con- 
trols. Taken together these characteristics (Table 6) merit fur- 
ther investigation of WARM for clinical use with [ n C]PIB, 
and the method may be equally applicable to other washout 
tracers. 

THEORY 

This section describes theoretical considerations regarding the 
three different approaches to the determination of the binding 
potential of [ n C]PIB (BP ND ) or DVR. 

REFERENCE REGION RATIO MEASURE (SUVR) 

The [ n C]PIB retention can be calculated by determining the 
accumulation relative to a reference tissue to obtain a ratio 
measure (SUVR). The ratio measure is the fraction of the region- 
of-interest integral of [ n C]PIB accumulation at steady-state, 
assumed to have been established no later than this time after 
injection (f s = 40 min), extended to the end (t e = 60 min). The 
volume of distribution of the tracer, Vj and the volume of dis- 
tribution of non-displaceable tracer (Vnd) define the binding 
potential BPnd- 

Vt = V ND + BP ND V ND (8) 

BPnd = r^- - 1 (9) 
Vnd 

The signal measured by PET is the accumulated mass of the tracer 
(Vt) relative to the integral of c a (t), the radioactivity of the tracer 
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FIGURE 7 1 The figure shows the average washout index © registered average MRI brain, the cyan insert shows the similar 

image calculated over the time interval T= 1890s for all reference region for the way compound as found in Moller et al. 

subjects. Mean values >1660s, superimposed on a non-linear (2009). 



Table 6 | Characteristic feature of the three methods used. 





WARM 


SRTM 


SUVR 


Method approach 


Calculation 


Non-linear 


Calculation 






regression 




Simplicity 


Moderate 


Complex 


Simple 


R1 correction 


Yes 


Yes 


No 


Noise-reducing 


Integration 


Fitting 


Steady-state average 


Reported parameters 


BPnd 


BPnd 


DVR 












" CBF" 


k2 





in the arterial circulation i.e., 



Vr 



ftm(t)dt 
ftc a (t)dt 



(10) 



The variable of major interest to [ C]PIB binding to A{5 for- 
mations is the quantity of specific (i.e., displaceable) binding. This 
binding is expressed as the binding potential (5Pnd)> the dis- 
placeable tracer relative to the amount of non-displaceable tracer 
in the tissue. In any given region of interest (ROI), an area of 
reference may be referred to as a region where equal amounts of 
non-displaceable tracer enter and leave the tissue. Such a reference 
fulfills the two requirements that, 

(1) tracer enters the tissue in proportion to blood flow, 

(2) unbound tracer (i.e., tracer dissolved in Vnd) clears the tis- 
sue at the rate determined by the blood flow, as defined for 
a tracer subject to flow-limited exchange across the blood- 
brain barrier. 



For the reference region, Vj = Vnd> an d Equation (10) defines 
the steady-state volume of distribution, 



ftm ND (t) dt 

ft C a (!) At 



(ID 



As the total amount of radioactivity in the arterial blood are the 
same in the two regions, i.e., the ROI (Equation 10), and the refer- 
ence region (Equation 11), the binding potential is obtained from 
the Equation (9). 



f 'm(t) dt 
BPnd = - 1 



(12) 



Equation (12) is valid when requirements (1) and (2) are ful- 
filled. However, for a tracer that disappears rapidly from the blood 
stream, (1) is invalid after the initial wave of tracer has passed 
the tissue. When the concentrations c a {t) in the circulation pri- 
marily depend on the wash-out and hence mainly on the CBF. 
Requirement (2) depends on the regional CBF and is valid for 
tracers subject to flow-limited exchange with brain tissue. 

It is improbable that we would find a reference region 
for all tissues and tracers that uphold these requirements, 
unless additional tomography is completed with a similar tracer 
of no specific binding and hence no displacement. However, 
when we identify a single homogeneous reference region for 
[ n C]PIB, we can mimic the tracer's behavior as reference for any 
given ROI. 

INITIAL [ 11 C]PIB DISTRIBUTION AS SURROGATE CBF MEASURE 

The two factors that affect the accumulation of [ C]PIB after 
the initial distribution are the magnitude of blood flow, which 
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mediates washout, and the degree of binding, which limits the 
rate of washout. In order to distinguish these factors and hence 
to assess the binding of [ n C]PIB in any given region or voxel, we 
determine the relative flow ratios R\ as the magnitude of CBF in 
a region or voxel, relative to CBF in the reference region, defined 
asRi =CBF/CBF ND . 

Cerebral blood flow can be assessed in two different ways, 
either directly by means of PET with [ ls O]water (CBF), or a 
surrogate measure (CBF PIB ) approximated indirectly from the 
[ n C]PIB signal at the peak of distribution. This approximation 
is made by first excluding the instantaneous arterial distribu- 
tion in the first 2 min after injection. In the subsequent 8 min, 
the initial distribution of the [ n C]PIB signal depends largely on 
the wash-in. Until maximum peak values are reached at At p , 
the wash-in can be expressed as the unidirectional clearance 
{K\) from the blood into the brain when no tracer has left the 
brain yet. 

We recently determined the permeability-surface product (PS) 
for [ U C]PIB (Gjedde et al, 2013). There the Ki measure from the 
initial distribution of the [ U C]PIB signal is related to the CBF by 
the Renkin-Crone formula (Crone, 1963; Renkin, 1964), 



= CBF^l - e~cm^ 



(13) 



where by definition the extraction fraction E from blood into 
tissue is, 

E=(l-e~^\ (14) 
From Equation (13), it follows that the wash-in ratio 



CBF PIB from the [ U C]PIB images. This derivation requires cal- 
ibration to a normal material of estimates of cortex Ki and 
CBF values, with tight adherence to the protocol used in the 
calibration. 

FLOW DEPENDENCE OF SPECIFIC BINDING MEASURE 

With the regional flow ratio measure Ri derived from [ ls O]water 
or [ n C]PIB analysis for each voxel with a signal n(t), we simulated 
flow-adjusted reference curves (mndW) for the corresponding 
voxel of the image, determining the tracer washout by the actual 
flow measured in the voxel. Here, the term mnd(0 represents 
the simulated dynamic wash-out of the tracer that would have 
been recorded from an individual voxel in the absence of bind- 
ing. The result is a simulated image of the dynamic passages of 
the tracer through every voxel as functions of time in the absence 
of any binding in any voxel. The simulation is based on the 
dynamic behavior in a fixed reference region (cerebellar gray mat- 
ter) given by mnd(0> which is the dynamic time-activity curve of 
the reference region. 

In order to determine the simulated reference curve «nd(0 
from the measured reference curve m^£)(t) , let a non-binding 
voxel or region have a flow CBF„, and let it refer to a measured 
region with flow CBFnd- Each region has washout rates given by 
fc 2 = CBF„/Vnd and fc^ D = CBFnd/ Vnd> respectively. Then _Ri 
relates to the variables as indicated by the following equations, 



H(0) 



CBF„ 



m ND (0) tf D CBFnd 



(17) 



where each voxel is the site of mono-exponential washout given 
by rate constants of the magnitudes, 



R 



PIB 



Ki/K[ 



ND 



(15) 



relative to a reference region is linked to the similar ratio for blood 
flow_Ri = CBF/CBFnd by the relative extraction fraction relative 
to the extraction fraction of the reference region £nd 



-.pib 



-End 



(16) 



ref 
P ' 

A tentative absolute measure for K\ can be estimated by 
additionally accounting for differences in the delay of the 
maximum peak and by normalizing for weight and dose, 
(m(Afp) /Atp) * Atp ei * weight/dose provides an delay nor- 
malized estimate of the magnitude of K\ . The estimate of K\ is 
scaled by a site- and tomograph-specific constant K s [ te for cerebral 
cortex values in order to derive tentative absolute flow estimates 



t \m ND (0) 



1 , ( «nd(0 
In 



This relation indicates that the R™ approximates R\ for 
extraction fractions similar to the reference region, which is the 
case for compounds with sufficiently high PS products and for 
regions where the CBF is close to the CBFnd of the reference 
region. With these limitations, and for want of known CBF val- 
ues, R\ lB may be estimated from the maximum signal intensity 

value m(Atp), such that _R PIB = m (Atp) /otnd ( Af p ef )> relative 
to the maximum signal intensity of the reference region at time wnere 



and 



t \"nd(0) 
where R\ relates the hi terms as, 

* /mNDW\ 
t \m ND (0)J 

such that Equations (19) and (20) yield «nd(0 as, 



In 



«ndW 
«nd(0) 



R, In 



«nd(0 = «ND(0)e 1 



WNPft) 

m N D(0) 



Und(°V 



or, 



«ND(f) = «ND(0) 



wnd(0) 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 
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where mnd(0 is the simulated flow- adjusted reference curve cor- 
responding to the measured curve n(t) when equal amounts 
of tracer enter, i.e., «(0) = «nd(0). When expressed for every 
voxel, the procedure yields a new dynamic reference image «nd (0 
that enables the comparisons among the flow-involved SUVR 
method and the flow-devolved WARM and SRTM methods, as 
the flow-devolved version of Equation (12) can now be expressed 
as, 



the [ n C]PIB signal depends only on the washout rate. The total 
volume of distribution in the binding region then is simplified to 

(Lammertsma and Hume, 1996; Tvfoller et al., 2009) 



where 



Vt 



k 2a k 2 



(29) 



f e n(t) dt 

BPnd = - 1 

f e n ND (t)dt 



(24) 



Figure 2 shows an example of the real and simulated non- 
displaceable (unbound) tracer time-activity curves for a small 
white matter region and a small putamen region, as well as the real 
time-activity curves for the ROIs. The variation in the simulated 
flow corrected reference curves, illustrate that a single uncor- 
rected reference curve may bias the result significantly and may 
only be valid for regions where the flow is equal to in this case the 
cerebellum. 

WASHOUT ALLOMETRIC REFERENCE METHOD (WARM) 

In the case of negligible input from the circulation after the initial 
brief uptake, the differential equations (25) and (26) 



and 



dm* (t) 
dt 



dt 



Kic a (t)-k 2a m*(t) 



K? D c a (t) 



ND 



(0 



(25) 



(26) 



are linked only while the tracer is dispersed from well-defined c a , 
i.e., during the brief uptake period until maximum peak (within 
2-10 min timeframe) when washout is assumed to be negligible. 

The term K\ is the unidirectional clearance of the tracer c a 
by the tissue, K^ D is the clearance of the tracer c a by the ref- 
erence region, m* and mj^ D are the measured PET signal in the 
tissue (with displaceable binding) and reference, respectively. The 
term k 2a defines the apparent measurable washout rate constant 
for the ROI. The term k 2 is the unknown washout rate for non- 
specifically bound tracer of the same region of interest, and fc^ D 
defines the measurable washout rate of non-specifically bound 
tracer in the reference tissue into the plasma. The uncoupling of 
the first and the second term on the right hand side of the equa- 
tions means that elimination of the first K\ and Kf D terms yields 
the equations. 



dm*(t) 
dt 

dm^ D (t) 
dt 



-k 2a m*(t) 



(27) 
(28) 



This defines the first-order decay for 1 -compartment first order 
kinetics, assuming that association and dissociation is sufficiently 
rapid in the tissue compartments. The Equations (27) and (28) 
each predict a mono-exponential washout from the time to where 



k 2a 



(1+5P ND ) 



(30) 



As the differential Equations (25) and (26) are linked in the initial 
phase (described by Ki ), the ratio _Ri = Ki /Kf D accounts for the 
difference of delivery to the regions of interest and reference. For 
the non-specifically bound tracer, we assume that the volume of 
distribution is the same in all regions, i.e., 



Ki 
k 2 



K 



ND 



(31) 



Consequently the rate constants k 2 and k^ D of washout of non- 
specifically bound tracer from tissue to plasma are similarly linked 
byRi, 

k 2 = R^ D (32) 
by combining which Equations (30) and (32) BP^d is found as 



BPnd = Ri 



k 2a 



1 



(33) 



The late uncoupled mono-exponential differential equations (27) 
and (28) can be rearranged to 



dm* it) 
m*(t) 

dm* UD (t) 

Integrating both sides the gives 
' m*{t) 



-kindt 



-kTdt 



In 



m*(0) 



In 



'ND 



(0 



'ND 



(0) 



-k 2a t 



-k^t 



(34) 



(35) 



(36) 



(37) 



where m*(0) = m*(to) and m^ D 
value at start of the washout to = 
(37) by (36) 



(0) = m^ D (fo) are the initial 
= 2 min. By dividing Equation 



In 



/ m NpM A 

, ( m*(t) \ 
l \m*(0)) 



iND, 



k 2a t 



(38) 
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The right hand side of this equation can be expressed in terms of 
the binding potential by use of Equation (33). 



In 



(m*(t)\ 

1 \m*(0)J 



(1 + BP ND ) 



(39) 



where (1+BPnd) is the distribution volume ratio(DVR). 
Alternatively written 



£i In 



BP ND + 1 



In 



V m ND(°)/ _ V ' V m ND(°)/ . 



(40) 



where (1 + BPnd) is the distribution volume ratio(DVR). 

Equation (39) can also be expressed as an allometric relation- 
ship between the logarithmic of the fraction of remaining f tracer 

in a ROI and the reference region where (^p) is the scaling 
exponent, found in this formulation by log-log linearization as, 



In 



'ND 



it) 



'ND 



(0) 



/DVR\ /m*(t) 
V^rTJ n \m*(0) 



(41) 



For direct calculation without linearization, Equation (39) was 
expressed in terms of integrations and differences of the loga- 
rithms, assuming constant BPnd and R\, which greatly reduces 
noise in BPnd estimates. 



BP ND (D = ^/o>K D (0)-lnK D (0)) * _ l (42) 
f Q (In (m*(f)) - In (m*(0))) dt 

That we further simplified with Equation (17) to the operational 
equation, 



BBnd(T) 



m*(0) / Q r (ln (mfojt)) - In (m^g)))) dt 
Kjd(°) /o T ( ln ( m *( f )) - In (m*(0))) dt 



1(43) 



When log transformed, the fraction (i.e., DVR) part of this equa- 
tions states that the nominator is the accumulated log-signal for 
the reference tissue relative to how much was present before 
washout, this difference is scaled by the initial tracer amount of 
the ROI. The denominator describes the accumulated log-signal 
for a ROI or voxel relative to how much was present before 
washout. This difference is scaled by the start amount of the 
reference region. Thus the fraction is corrected both for flow, 
i.e., initially deposited tracer, and the exponential behavior of 
the washout. Figure 1 illustrates the behavior of the nominator, 
denominator and BBnd(T) of Equations (43) and (40) for simu- 
lated ROI and reference curves with 20% added Gaussian noise. 
As seen the nominator and denominator from Equation (40) (in 
panel 2 from the left) are stabilized by the integration in Equation 
(43) (panel 3 from the left). Panel 4 illustrates how the BB ND (r) 
values converge towards the theoretical result. 



SIMPLIFIED REFERENCE TISSUE METHOD 

The simplified reference tissue method (SRTM) (Lammertsma 
and Hume, 1996) yields binding potential when a single tis- 
sue compartment model fits the data. SRTM solves differential 
equations similar to (25) and (26) (equations (1) and (5) in the 
paper of Lammertsma and Hume (1996)). The method assumes 
that these differential equations are coupled by tracer distribu- 
tion in the circulation throughout the entire duration of the 
tomography. 
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